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The covariant gauge invariant perturbation theory of scalar cosmological perturbations is devel- 
oped for a general Scalar- Tensor Friedmann-Lemaitre-Robertson- Walker cosmology in a vacuum. 
The perturbation equations are then solved exactly in the long wavelength limit for a specific cou- 
pling, potential and background. Differences with the minimally coupled case are briefly discussed. 

I. INTRODUCTION 

In the last few years a tremendous advance in the accuracy of cosmological observations has revolutionised our 
understanding of the dynamics of the Universe. In particular, it has been realized that the universe is evolving in a 
way that is incompatible with what one would expect in a homogeneous and isotropic baryon dominated universe. One 
way of dealing with this problem is to suppose that the Universe is still homogeneous and isotropic but dominated by 
a different form of energy density: the so called Dark Energy (DE). Many different ideas regarding the nature of DE 
have been proposed in the last few years, ranging from the cosmological constant [l| to Quintessence Chaplygin 
gas and Phantom energy [3| , but so far none of them has provided a fully satisfactory solution to the problem. 

Another interesting approach is based on the hypothesis that on cosmological scales gravity works in a slightly 
different way than in General Relativity (GR). In this way DE acquires a geometrical character, i.e. some of the 
observations can be explained by having a different behaviour of gravity on cosmological scales, rather than a new 
(and as yet undiscovered) form of energy density. 

The idea that gravity might work differently on different scales is suggested by many fundamental schemes such as 
quantum field theory in curved spacetime Q, as well as the compactification of internal spaces in multidimensional 
gravity [l,0|i the low energy limit of string theory and some braneworld models [9]. Furthermore, GR has only been 
well tested on small scales and at low energies, so currently we lack strong constraints in the other regimes. In this 
respect cosmology is an ideal testing ground because it offers an indirect way to test GR on scales and at energies which 
are not achievable in Earth-bound laboratories, but were realized throughout the history of the universe. Therefore, 
investigations of cosmology in alternative theories of GR has a double relevance; it contributes to the understanding 
of DE and provides a testing ground for GR on different scales. 

A key step in investigating the cosmology of an alternative theory of gravity is the development of a full theory of 
cosmological perturbations. This will make it possible to use some of the best and most precise sets of observational 
data currently available, such as the WMAP observations of the Cosmic Microwave Background (CMB) anisotropics, 
allowing us to understand the viability of these theories and constrain their parameters. 

In this paper, we will develop the theory of linear scalar perturbations for one of the most studied extensions of 
GR: Scalar- Tensor theories of gravity. This will be done using the 1-1-3 covariant and gauge invariant approach, 
developed in [l3| and successfully applied to GR in the presence of minimally coupled scalars field [ll|. Although 
this formalism can be applied to perturbations of any spacetime, our work will be focused only on almost Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universes which correspond to standard linear perturbations of FLRW models 
(see for details) and in absence of standard matter. The non-vacuum case will be treated in a forthcoming paper. 
Scalar field dominated universes in the Scalar- Tensor framework have attained prominence through the Extended and 
Hyper- extended InHsLtion scenarios [Tsl. fl^. so an analysis of the evolution of perturbations is important for obtaining 
a complete description of structure-formation in these models. In addition, the general formalism presented here could 
be used in situations different from inflation in which a scalar field dominates. As in [Tl| . emphasis is given here to 
curvature perturbations [Tsj . which are naturally gauge invariant, rather than metric perturbations p^ . |l3,[ll which 
play no explicit role. Bardeen's formalism has been applied to this situation in specific cases of coupling, potential 
and background [l9| . There has also been quite an extensive analysis of perturbations in non-standard theories of 
gravity given by Hwang in (20| . The aim of this paper is to give the general perturbation equations for a general 
coupling, potential and (FLRW) background using the full covariant approach to cosmologyand choosing a set of 
variables which are as much as possible of straightforward physical interpretation. 

The paper is organised as follows. In section |TT1 after a few preliminary remarks on Scalar- Tensor gravity, we set 
up the formalism based on the natural slicing of the problem {(jj = const} and on its geometric characterisation 
through the unit vector which is orthogonal to these surfaces. We also characterise the effective thermodynamics 
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of the non-minimally coupled field by treating it as an effective fluid. In section IIIII we specify the general features 
of the background and define the key gauge invariant dimensionless variables. In section ITVl we present the general 
perturbation equations, and, after specifying the coupling, the potential and the background, we give an exact solution 
of the perturbation equations in the large wavelength limit. In section 5, we work out a specific example giving an 
exact solution in the long wavelength limit i.e. on super-horizon scales. 

Unless otherwise specified, natural units {h ~ c = ks = SnG = 1) will be used throughout this paper, Latin indices 
run from to 3. The symbol V represents the usual covariant derivative and d corresponds to partial differentiation. 
We use the —,+,+, + signature and the Riemann tensor is defined by 



where the W§j^ is the Christoffel symbol (i.e. symmetric in the lower indices), defined by 

1 



r.?"*^ {gbe,d + ged.b — gbd,e) 



(1) 



(2) 



The Ricci tensor is obtained by contracting the first and the third indices 

Rab — g'^'^Racbd ■ 

Finally the Hilbert-Einstein action in presence of matter is defined by 



(3) 



A 



(4) 



II. PRELIMINARIES 



A. Scalar Tensor Theories of Gravity 



Scalar- Tensor theories of gravity involve a new degree of freedom of the gravitational interaction in a form of a 
scalar field non-minimally coupled to the geometry. These theories were first proposed as a modification of GR 
able to completely include Mach's principle From a technical standpoint, this can be achieved by allowing the 
Gravitational constant G to vary throughout the spacetime, and in particular to be a function of a scalar field (so 
that frame effects are avoided). However, such a scalar cannot be found within GR because they can contain first 
derivatives of the metric and fall off more rapidly than r~^ . As a consequence, G is imagined to depend on an 
additional scalar field 0. 

Like many of the extensions of GR, scalar-tensor theories appear in many fundamental schemes and have been 
proposed as models for dark energy because their cosmol ogy naturally leads to the phenomenon of cosmic acceleration, 
which is the characteristic footprint of dark energy [l^, [2^ . 

The most general action for Scalar Tensor Theories of gravity is given by (conventions as in Wald [13]) 



I dx 



-g 



\F{<i>)R - iVa0V"0 ~ y(0) + 



(5) 



where V{(j)) is a general (effective) potential expressing the self interaction of the scalar field and Cm represents the 
matter contribution. 

Varying the action with respect to the metric gives the gravitational field equations: 



F{<jy)Gab = F{<jy) ( Rab 

"^aCf^^bCt) 



T. 



ab 



gab 



gabR 



(VbVa - gabVcV') F{<j>) , 



while the variation with respect to the field (j) gives the curved spacetime version of the Klein - Gordon equation 

V,V>-HiF'(0)i?-l/'(0) = O, 



(6) 



(7) 



where the prime indicates a derivative with respect to 0. Both these equations reduce to the standard equations for 
GR and a minimally coupled scalar field when -F((^) = 1. 
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Equation ([6]) can be recast as 



where Tj'j has the form 



1 



Gab 



1 



F{<P) 



T. 



ah 



Va(bVb(b - gab ( 2 Vc0V^0 + 1/(0) ) + VbV„F(0) - 5afcVeV^F(0) 



Provided that (j)^a 7^ 0, equation ([7]) also foUows from the conservation equations 

^'T^b = . 



(8) 



(9) 



(10) 



The fact that the gravitational field equation can be written in this form is crucial for our purposes. In fact, the 
form of ([H]) allows us to treat scalar tensor gravity as standard Einstein gravity in presence of two effective fluids 

with energy momentum tensor ■j^^ and T^j^. This implies that, once the effective thermodynamics of these fluids has 
been studied, we can apply the standard covariant and gauge-invariant approach to cosmology. In this paper we will 
consider only the vacuum = case, leaving the non- vacuum case to a future paper. 



B. Kinematical quantities 



In order to give a proper description of the kinematics of the effective fluid associated with the scalar field we have 
to assign a 4- velocity vector to the scalar field itself. Following 11, 25j we assume that the momentum density 
is timelike: 

Va0VV<O, (11) 

in the open region A of spacetime we are dealing with. This has two important consequences: Firstly (/) — const. 
specifies well-defined surfaces in spacetime and secondly, these surfaces are spacelike (as consequence of (dU). We can 
therefore choose the 4-velocity u° to be the normalised timelike vector 

y- = _Z!^^ U-Ua^-l, i> = ^^u''Va<jy={~^a<jy^''<l>Y'\ (12) 

■0 

where -0 = denotes the magnitude of the momentum density (simply momentum from now on). The vector (jl2p is 
timelike because it is parallel to the normals of (i.e. orthogonal to) the surfaces {0 const.} and it is also unique 
because, since these surfaces are well defined, their normals are unique. 

Note that this definition might not appear well posed in general because could oscillate and change sign even in 
an expanding universe. However, equation (|lip and the regularity of the behaviour of (see |ll| for more details) 
ensures the possibility of an extension by continuity. 

Once the velocity field has been chosen, the projection tensor into the tangent 3-spaces orthogonal to the flow 
vector can be defined as: 

hab = gab + UaUb ^ h'' bh^ ^ = ^"c , KbU^ = , (13) 

and the derivation of the kinematical quantities can be obtained by splitting the covariant derivative of Ua into its 
irreducible parts p^ : 

VfcUa = VfeUa - ^aUfc , ^ bUa ^ ^Qhab + <Jab + ^^ab , (14) 

where is the totally projected covariant derivative operator orthogonal to u", = Ua is the acceleration [AjjU^ — 
0), Q is the expansion parameter, aab the shear (cr^a = crabU^ = 0) and LOab is the vorticity {ujab ~ ^^[afc]; ^abu'' = 0). 
Important consequences of the choice ([T2]) are 



Va0 = 0, 



(15) 
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i.e. the surfaces {(j> = const.} are the 3-spaces orthogonal to and 

e = -v<,(^) = -i [F'(<^)i?-v^'(^) + V'] ^0 

(Tab = -^ha^W^V(, [Vd)0] + i/lafcVc(^) , 



(16) 
(17) 
(18) 



where the last equality in follows on using the Klein - Gordon equation ^ . We can see from that V' is an 
acceleration potential for the fluid flow [l3|. Note also that, as consequence of (fT5)) . the vorticity vanishes: 



= 



(19) 



which implies that Va is the actual covariant derivative operator in the 3-spaces (j) = const . Finally, it is useful to 
introduce a scale length a along each flow -line so that 



a 6 



(20) 



where H is the usual Hubble parameter if the Universe is homogeneous and isotropic. 



C. The non-minimally coupled scalar field as an effective imperfect fluid 



Using the definition of the 4-velocity given above, we can deduce the thermodynamical properties of the effective 
fluid described by the non minimally coupled scalar field. The general form of the energy momentum tensor ([9|) is: 



where the energy density /i^ and pressure p^ of the scalar field "fluid" are: 

1 



IJ'4> 



F(0) [2 



+ V{^) - QF{4>) 



(21) 



(22) 



P0 



1 



F{4>) 



1 



and the energy flux and anisotropic pressure are given by: 



la 



Aa 



(23) 



(24) 



<7ab ■ 



(25) 



Since the anisotropic pressure iVab 34^ and the heat flux qa are different from zero, the effective fluid associated with 
a non-minimally coupled scalar field is, in general, imperfect. The form of (j24p and (|23p shows explicitly that this is 
a direct consequence of the non-minimal coupling and that, when </> is minimally coupled (i.e. F{(l)) — 1 — constant) 
we have qa = = i^ab- This is in agreement with the results of ll[, in which a minimally coupled scalar field is 
considered to be a perfect fluid. It is also worth noting that, since iTab and qa depend on the acceleration and the 
shear, when we consider FLRW universes, where these quantities vanish, the effective fluid behaves like a perfect fluid 
[3^ . However this is not true in the perturbed universe. 

The twice contracted Bianchi Identities lead to the energy and momentum conservation equations: 



fi + e{fi+p)+ -Kabo"^ Vag" + aaq°- = 



(26) 
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and 

Aa(p+p)+^aP + h\(WhTr''c + qc) + (^^r'a + ^e/i^^ gfa = . (27) 
If we now substitute p and p from (|22p and ((23)l into ((26)l . we obtain, as expected, the Klein - Gordon equation ([7]): 

+ e0_iF'(0)i? + y'(<^)-o, (28) 

an exact ordinary differential equation for in any space - time with the choice (jl2p for the four - velocity. With the 
same substitution, (j27p becomes an identity for the acceleration potential ip. 

In our calculation the ratio between the pressure p and energy density n of the effective fluid defined above is given 
by 



pF{(b) 



^2 + ieF(0)+F(0) 



(29) 



This is a useful parameter of the theory, but should not be confused with the barotropic factor typical of perfect 
fluids. In the same way, the quantity will be denoted by 



fj- 

c — — + w , (30) 
w fi 

even if it does not represent the proper speed of sound. 

III. GAUGE -INVARIANT PERTURBATIONS AND THEIR DYNAMICS 

A. Background dynamics 

The flrst step in developing a theory for the evolution of cosmological perturbations is the definition of the back- 
ground. For this purpose it is useful to employ the 1+3 covariant approach to cosmology jK|. This method consists 
of writing the Einstein Field Equations and their integrability conditions as a system of six exact evolution and six 
constraint equations {1+3 equations) for a set of covariantly defined quantities (which include the expansion, the shear 
and the vorticity, already defined in section FlIAp . The advantage of such a re-parametrisation is that the treatment 
of the exact theory is considerably simplified, making it easier to find background solutions, even in more compli- 
cated cosmological models (like the Bianchi spacetimes). In addition, many of these 1-1-3 variables can be shown to 
be gauge-invariant [27^ and it is these quantities that form the building blocks of our gauge-invariant perturbation 
theory. 

In this paper we limit our analysis to the study of perturbations of spatially homogeneous and isotropic cosmological 
models (FLRW models), in which only a spatially homogeneous non-minimally coupled classical scalar field is present. 
This means that in the background we have 

a = w = 0, Va/ = 0, (31) 

where / is any scalar quantity; in particular 

VaM = V„p = => V,V^ = , A, = . (32) 
In this way the 1+3 equations reduce to a system of only two equations: the Raychaudhuri equation 

e + ie2 + i(/i^-h3p^) = 0, (33) 



and the Gauss - Codazzi equation 



R=2 



(34) 

where R represents the 3-ricci scalar of the (j) = const surfaces, and is closed by the conservation equation (j26p . 
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Expressing ((26)) . ((28|) . ([341) . (|33|) in terms of the Hubble parameter and the momentum variable ip, we obtain the 
background (zero -order) equations: 



3H + = 



1 



3ij2 + SiC = 



(35) 
(36) 



+ 3Hil; - ^F{(P)R + V'{(j)) = ^ (37) 

A + 3i/ (^v^2 + ^ei^(0) + i^((^)^ , 

where all variables are functions of cosmic time t only, and 

R = QK/a^, K = const (38) 
as required by the properties of homogeneous and isotropic spacetimes. 

B. Gauge - invariant perturbation variables 

1. Spatial gradients 

Using the 1+3 variables we can define exact quantities that characterise inhomogeneities in any space -time, and 
also derive exact non- linear equations for them. For example, to characterise the energy density perturbation, the 
natural variable is Xa — V a^J'■ This vector represents any spatial variation of the energy density /i (i.e. any over- 
density or void) and it is in principle directly measurable [28| . However, a more suitable quantity to describe density 
perturbation is 

Va = -Xa - -VaM , (39) 

where the ratio Xa/ allows one to evaluate the magnitude of density perturbations relative to the background energy 
density and the presence of the scale factor a guarantees that it is dimensionless and comoving in character. The 

magnitude of T>a is closely related with the more familiar quantity — , the difference being that T> ~ (VaV^y/^ 

M 

represents a real spatial fluctuation instead of a fictitious gauge dependent one. In fact, it is easy to show that the 

S ^ 

Bardeen variable em which corresponds to — in the comoving gauge, is the scalar harmonic component of Va (i.e. 

its scalar "potential"). 

Other important quantities relevant to the evolution of density perturbations are 

= aV,e, Ca = a^VaR, (40) 

which represent, the spatial gradient of the expansion and the spatial gradient of the 3-Ricci scalar respectively. 
These vectors, together with Va, vanish in the background FLRW model and are therefore gauge-invariant by the 
Stewart- Walker lemma p9j . 

At this point one could derive the exact non - linear evolution equations for density perturbations from the definitions 
above and the equations ([26]) . (j28|) . ([34]) . (I33|l . However, since in our case we are considering an effective fluid 
representing the non-minimally coupled scalar field, we may want to characterise directly the inhomogeneity of 4>. 
This cannot be done using Waip because our frame choice implies that the spatia,l gradient Vafj) identically 
vanishes in any space -time. As a consequence (and just like the minimally coupled case pU. [30j). the natural gauge 
invariant perturbation variable for the inhomogeneities of the scalar field in our approach is the spatial variation of 
the momentum Vqi/'- Following the same reasoning used for Va we can define the dimensionless gradient 

^a = T^aV- • (41) 
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Since an exact FLRW Universe is characterised by the conditions (|3ip . ((32)) . vanishes identically in the background 
and therefore is gauge - invariant. In addition, by comparing (|18p and (j4ip . we see that is proportional to the 
acceleration (ai, = '^b/a^). Thus, this quantity can be interpreted as a gau ge in variant measure of the spatial variation 
of proper time along the flow lines of u° between two surfaces (j) = cons^. [30l[3l|. 



2. Linearization 



Since the variables (T>a, Za, Ca, ^'a) are completely general, solving their evolution equations in full would be 
equivalent to solving the full Einstein equations, a very difficult task. On the other hand, the current observations 
suggests a universe which is very close to FLRW. This means that we can limit our analysis to situations where the 
real universe does not differ too much from the background, i.e. when the magnitude of the quantities in (PTjl . (|32p . 
([39)1 and (|40|) are small. This can be achieved by treating the quantities that do not vanish in the background as zero 
order and those that vanish (and are hence gauge-invariant) as first order, retaining only terms which are first order 
in the gauge-invariant perturbation variables. 

In our specific case (background given by FLRW universe -I- a scalar field), the zero order quantities are given by 
Q, /i, p and the first order ones are aab, QajT^ab, T^a, Za, Ca, ^'a- So for example, when we linearise the equations, we 
drop the vector aabD'^ because it is second order, but we keep the quantity OVa, where it is understood that O is the 
expansion in the background. 

In the linear approximation we find two relations between the perturbation variables defined above: 



Ca = -^&a^Za + 2fia^Va 



obtained by taking the gradient of ([M)) and 



2F(0) - 3F(,/))] Va = ^ e/^(^)^ _ 3F(</))Ca , 



obtained from the definition of the effective energy density. Using these equations we can express the heat flux qa and 
the anisotropic pressure TTat in terms of these variables. For the heat flux we have 

F(0) F(0) ^a 



F{(j)) F{(l)) a2 ' 

where we have also used the relation between acceleration and "^a- For the gradient of the anisotropic pressure we 
have 

-. , ^FW ( _C^ _ M^i _ EMllA 
F((^) \2a^Q a^e F(0) j ' 

where we have also used the linearised vorticity free shear constraint 

- aVVab - ag, = . (42) 



3. Scalar gauge - invariant variables 



The vectors {'Da, Za, Ca, ^'a) contain information about the evolution of the energy density perturbations that is 
not necessarily related with the formation of local inhomogeneities [l2| . The relevant parts can be extracted using a 
local decomposition. For example in the case of the gradient of Va we have: 

aVfePa = Aab = ^habA + Y.ab + Wab , (43) 

where 

= A(^f,) - ^A/lab , Wab = A[ah] , (44) 

are the symmetric, trace - free and anti - symmetric parts of Aaf, respectively. The first represents spatial variations of 
2?a associated with change in the spatial anisotropy pattern of this gradient field and it can be related to the formation 
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of non-spfierically symmetric structures. Tfie second one represents spatial variations of T>a due to rotation of the 
density gradient field. Using the definition of Wab it is easy to see that this tensor is proportional to the vorticity 
[T^ and therefore vanishes for our choice of 4- velocity u"' . The first term A in (j43|) is the trace of l^ab- It represents 
the spherically symmetric spatial variation of the energy density /i and is the variable most closely associated with 
matter clumping. The evolution of this variable will be the focus of this paper. 

In the same way one can obtain the corresponding decomposition of the variables ()40|) and consider only the trace 
of their gradient. In this way we have the set of scalar gauge invariant- variables: 

A = aV"X>a , Z = oS/^-Za , C = aSI^Ca , * ^ oSI°"^a , (45) 

respectively giving the energy density, expansion, 3-curvature scalar perturbations and the spatial distribution of the 
momentum ^. Constraint equations analogous to (|42l) and ([47|) then relates Z, C, A and 4": 

C = -^Qa^Z + 2fia^A , (46) 



2&F{(j)) - 3F(0) A = 2a^Q + ©^(0)) * - iF{(P)C . (47) 

In the minimally coupled case {F{(j)) = constant) this last equation reduces to 

A = (u; + 1)^- , (48) 

which is equivalent to the result found in ^lH, and we conclude that, in this case, A characterises scalar field clumping. 
However, this is not true in general for a non-niinimally coupled scalar field because (j47p is related to the variable C, 
the details of the coupling function F and its derivative, as well as specific features of the background like 8 and the 
sign of (/). 

This is not surprising: A represents the spherically symmetric clumping of the effective fluid described by the stress 
energy tensor ^ and, in general, it has no direct relation with the inhomogeneities of the scalar field. However, the 
behaviour of this quantity can still be interesting. For example, in the context of the inflationary scenarios, even if the 
perturbations of the scalar field are described by 5*, only A is connected with the generation of the seeds of structure 
formation. 

Finally, for later convenience, it is useful to define the vector 

= ^^a^Va (49) 
and its comoving divergence $ defined in [TTj. This quantity is equal, up to a constant, to the combination 

$jv = ^ {^A ~ ^h) (50) 

of the Bardeen potentials ^a and and represents the Newtonian gravitational potential. We will consider the 
evolution of <f>Ar in the specific example given in the last section. 



C. Entropy perturbations 

As we have seen, even if the effective fluid in our background takes on the form of a perfect fluid, this is not the 
case in the perturbed universe. Furthermore, its equation of state does not have the simple barotropic form p = p{^), 
so that the perturbations are in general non-adiabatic and we need to introduce another variable to characterise the 
spatial variation of the entropy density s. The form of this variable is obtained by taking the spatial gradient of the 
equation of state p — p(/i, s): 

aVaS, (51) 

where we have the usual thermodynamic partial derivatives at constant density and entropy [SG*]. The second term 
in the equation above represents the comoving gradient of the entropy s: 

VaS , (52) 



aVaP = — 



~ dp 

OS 



. _ a dp 

■a ~ r^ 

p OS 
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and corresponds to an effective entropy perturbation variable. Defining tfie comoving fractional pressure gradient 

Pa = — VaP, equation (|5ip can be written as 
P 



which means that 



F — —V — P 

W 



w w 



(53) 



(54) 



where we have used the definition of the effective barotropic factor w — p/fj,. Note that, £a is proportional to Va if 
w — const and it is zero in the case of a perfect fluid (c^ — w). The definition of £a is very close to the definition 
of the Bardeen variable 77 p^ . i.e., representing the difference between pressure and density perturbations when the 
perturbations are not adiabatic. 

As for Va2?6, we can decompose Va^b in such a way to isolate only the scalar part that is relevant for our purposes: 



£ = a' 



-A 



w w 

Substituting for W'^w from ((29|) . we obtain, after a lengthy calculation. 



(55) 



F 



F 



C 



6tp'^-eF + 3F~3ij^F" 



F w i-j, A^w F a? Q 3w 11 F 

F (e^ F -3nF -SeP + 18^jj^ + 9f) - 3F {3^^ - OF + 3F 



+ 



%wQ^iF^ (1 + w) 

This will be will be used later in the derivation of the perturbation equations. 



A 



(56) 



IV. PERTURBATION EQUATIONS 



A. First - order linear equations 



We are now ready to write down the linear evolution equations for the scalar gauge-invariant variables defined 
above. Since we have four perturbation variables A, 4", C [s^l and Z, we will have a system of four differential 
equations. Using the constraints and (|T7)l . it can be reduced to a system of only two equations. However, in the 
following we will consider equations for the variables A, C, i.e., we do not use the constraint (|47p. This choice is 
motivated for the sake of simplicity: the system for these three variables is much easier to manage than the general 
equations for two variables only. The choice of C as our third variable is, of course, purely arbitrary: one could, for 
example, couple the first order evolution equation for A and 5" to the equation for Z, but we choose C because it 
makes it easier to make a direct comparison with the minimally coupled scalar field and barotropic fluid cases. 

Starting from the definition of the gauge invariant variables given above and using the fact that at linear order 

aV aXiy = a {v aXbj holds, we obtain the following system of equations: 



A = 



3 (w 

4 



DC 



wO u(w 

2' ^ 



1)0-1 



A - -VVag'' + -95/ :, 



(57) 



C = 



where 



3(1 + w) 
4a4e 



+ ^ A 



3{l + w) 



~7 , ^ 



3a-. 1 y^^^'4]s-2a^^^^.<^^^3K(4-^'-^ 
3 (1 -f w) ^ V / V a"^ B fc) 



(58) 



(59) 
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together with an equation for ^' given by 

6Fe/i {l + w)-9iJ,F (1 



c 



80292 (^11,2 ^QF^ 

iP^ (292 -3 (1 + 3w) Ai) + 12 9V'V'' + 6f^9M(l + 3w) -692F 
6 9 - ^) 

F (2 92 + 3 (2 + 3 w) m) F (9 (1 + w) m - 2 92) - 6 9^^ 



(2F9 + 3f) 



2 

Fa2 
V'2 - 9F 



39 - 9F 



5 



402 (^2 _ 



A/i 



OF 



{1 + w) fi (^ijr 



OF 



wF 



(l + w) [i^^-OF 
3k 



V' + 5 ) 5 



3fc 



a 



(1 + w) (^V^ 



9F 



9fcF 



4a2 9^ 



(V2 



9F 



and the constraint ([T71) : 



29F - 3F 



40^9 h/>2-9F 1'-3FC 



(60) 



These equations reduce to those for a minimally coupled scalar field when F((j)) — 1. Substituting for the heat flux, 
the anisotropic pressure and the entropy in ([57|) . ([58]) and ((60|) we obtain: 



A = C + 6a A + Ca + I?A V^C + £a V^A + J^a V^^- , 
C= K{AcC + BcA + Cc'i>) + {Vc C + £aA + J^c ^') , 
* = C + 6^. A + * + P^/ V^C + £^ V^A + V^^- . 



(61) 



where for sake of simplicity the time dependent coefficients have been indicated with the curly letters A,... T. 
Expression for these coefficients are given in the appendix. Note that as in the minimally coupled case, the evolution 
of either the momentum perturbation ^ or the density perturbation A is closed at second order (in time derivatives). 



B. Harmonic decomposition 

((6T|) is a system of partial differential equations which is far too complicated to be solved directly. For this reason we 
follow the standard procedure and perform a harmonic decomposition of the perturbation equations This operation 
is usually performed using eigenfunctions of the Laplace-Beltrami operator on the 3-surfaces of constant curvature 
that represent the homogeneous spatial sections of the FLRW universes. 

The situation is different in our approach: we only consider quantities defined directly in the real perturbed universe 
and, in general, the spatial derivative Va is not a derivative on a 3-surface (unless the vorticity is zero [T^V Therefore, 
in order to use a harmonic decomposition we need to construct a new set of harmonics using operators constructed 
from these derivatives which are independent of proper time (i.e. constant on the fluid flow lines). This is easily done 
by writing the Laplace-Beltrami operator for Vq and defining our harmonics as the eigenfunction of the Helmoltz 
equation (see [l^l for more details): 

fc2 fc2 fc2 
V2g = -g, V^Qa- jQa, V^Qafc^ jQafa, (62) 

a2 

where fc = 2i:a/X is the wavenumber and Q = Qa = Qab = 0. 

Using these harmonics we can expand every first order quantity in the equations above, for example [38j . 



A^5]A(,)(t) 



(63) 
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where ^ stands for both a summation over a discrete index or an integration over a continuous one. Using the 
harmonics Q, together with equations ([5^ . the system can be written as 

Afc = Ma - + Sa - Afc + Ca - 



'J'fc = M* - Cfe + - Afc + - 



In this way we reduce the system to one of ordinary differential equations that describe the evolution of the inhomo- 
geneities for a given scale and can be solved in the standard way, once the background has been specified. 



1. The long wavelength limit \ ^ H 



Before calculating the evolution of perturbations for specific FLRW backgrounds it is worth considering an inter- 
esting sub-case. As we can see from the equations (j62p . by only considering modes with wavelengths much larger 
than the the Hubble radius, all the Laplacian terms can be neglected in (|6ip . From the second equation in (pT|) . this 
implies immediately that the curvature gradient Ca and its divergence C are conserved when the background universe 
is flat {K = 0). In this case the system (|6T]) can be simplified to give 

A= ^aIk=oC+6a|,,=o^+CaIk=o*' (64) 
*= ^*Ia-=oC+S*Ik=o A+C*|^^o*. (65) 

where C is constant in time. Note that in the above system, no spatial derivatives appear and as a consequence the 
spatial part of (A,C, ^') can be factored out leaving a second order system of ordinary differential equations. In the 
treatment of the evolution of perturbation on super-horizon scales and flat spatial geometry, we solve directly (j64p . 
using for simplicity a plane wave expansion for the spatial part of the perturbation variables. The long wavelength 
limit will be used explicitly in the example given in the next section. 



V. APPLICATION: THE CASE = -^(^^ l/(<^) = X<pP 

Let us now consider a specific example using a background which may be obtained, with a very natural choice of 
some integration constants, in the context of Noether symmetries applied to the FLRW metric In this model the 
coupling function is given by 

F(0) = -^(r)02, (66) 

and the potential 

^(0) = A(/)2f(''\ (67) 
where ^ and p depend on a single free parameter r, via 

' 12(r + l)(r + 2) ' ' 2, + 2r ^ ' 

Equation (|66p is different from the usual non-minimal coupling F{(f>) = 1-1- Ci^^ because in general it does not contain 
the Hilbert Einstein term. In addition, the value of r determines the sign of the effective gravitational constant 
Geff — F~^. Since the sign of the function ^(r) is not always the same, there are, in general, values of r for which the 
gravitational interaction is effectively repulsive. These facts deserve a detailed discussion that is outside the purpose 
of this paper (in which this model is just a worked example) and will be given in a forthcoming paper. However, 
roughly speaking, one can imagine that the real coupling is actually F{(j>) = 1 — ^0^ and that the condition is 
realized only at early times so that ordinary gravity is recovered at late times. 
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Using the Noether symmetry approach with the above couphng and potential leads to the solution 

a{t) = at" ; (f>{t) = /Jt"; K^Q, (69) 

with 

, , 2r2 + 9r + 6 , , 2r2 + 9r + 9 

nir] = — — ; m(r) ~ — — . (70 

^ ' r{r + 3) ^ ' r{r + 3) ^ ' 

which represents a specific case of a general exact analytic solution for the background equations [l^, [H, [s^ . De- 
pending on the value of r, (|69p can represent power law inflation behaviour (n(r) > 1) or a Friedmann-like phase 
(0 < n(r) < 1)[3^ and its character is essentially related to the choice of the form of the potential (via the choice of 
r). The parameter a is free while (3 is linked to A and r through 

(6 + r)(3 + 2r)2 

A — ^ — 3 3+2r 

8r(3 + r)2(2 + 3r + r2)'^ 

It turns out simpler to leave (3 free and derive A, but of course the real free parameter is this last one. Let us write 
down the time dependence of some important quantities. The expansion parameter is 

e(t) = 3^; 

the effective energy density is 

3(6 + 9r + 2r2)2 1 



r2(3 + r)2 i2' 

which is a non negative function of r regardless of the sign of Geff, and the effective barotropic index is 



2 2 

w — — - — 



3 6 + r(9 + 2r)' 

which is constant so that in this background the effective fluid behaves like a barotropic fluid [iot. 

Let us now write the perturbation equations in the long wavelength limit. As we have seen, in this limit the spatial 
dependence of the perturbation variable can be considered as a plane wave and can be factored out. Substituting the 
quantities above in the perturbation equations (j64p and in the constraint (j47p we obtain, after lengthy calculations, 

C 3(l + r)(4+r) 

A ^ nir)A+I{r)'^ + M{r)—t '-(3+.) , (71) 

a'' 

* = A/'(r)A + C'(r)1' +7'(r)-^r , (72) 
(7 = (73) 

C 3(l + r)(4+r) 

Q(r)A = Uir)^ +T{r)—t ^Ts+^ , (74) 
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where 



■H{r 

lir 

M{r 

N{r 

0{r 

V{r 
Q{r 
T{r 



54 + r (l53 + 117r + 23r2) 

r(3 + r)' ' 
-2 (3 + 2r) (6 + 5r) (6 + 15r + 4r2) 
r (3 + r) (6 + r (9 + 2r)) ' 
r (3 + r) (6 + 5 r) [21 + 2 r (21 + 12 r + 2r2 
6 (6 + r (9 + 2r))^ 
3 (6 + 9r + 2r2) (6 + 15r + 4r2) 
2r(3 + r)^ 
3 [12 + r (46 + 36r + 7r2)] 

r(3 + r)^t 
r [36 + r (90 + 57 r + 10 r^)] 
12 (6 + r (9 + 2r)) ' 
3 (6 + r (9 + 2 r)) , 

ir2(3 + r)' (3 + 2r) , 
2 (3 + r) (6 + 5 r) . 



Using the constramt ([Ti]) . the above system can be reduced to a particularly simple closed equation for 



(6 + r) 
" (3 + r) t 



r(3 + r) C 3(1+^) (4+r) 
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(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
(81) 
(82) 
(83) 

(84) 



where C is constant. Note that in the above equation the role of the spatial curvature perturbations C depends 
strictly on the value of r. In particular, the C term can grow or decrease in time leading to a different effect on the 
overall perturbation dynamics. 

Solving ([51)1 and using the constraint (|47p . we obtain the general solutions for 5* and A 



(3 + r)^ 
12 (12 + 6r + r2) ^2 



6 + 6 r + r-' 6 + r 



^ _ r2 (3 + r)^ (12 + r(21 + 2r(6 + r))) C ^_2^!i0^ 
18 (12 + r (6 + r)) (6 + r (9 + 2r)) ^ 
^ 2 (3 + r) (6 + 5 r) 6+. 

+ '^n — — -t ~ 

° 3 (6 + r (9 + 2r)) 



(85) 



(86) 



where '^q is an integration constant. 

The detailed properties of the physics associated with these solutions will be discussed in detail in a following paper. 
Here we will limit ourselves to making few brief comments. First of all in the above solution A has the same time 
dependence of ^ i.e. the scalar field clumps in the same way as the effective fluid. 

In the inflationary regime (r < — 3 — \/3, — 3<r<— 3 + -y/S, r > 0) the long wavelength perturbations are frozen 
and decay (see Figure [T|), unless —6 < r < —3 — V3 for which the second mode of ((85|l is growing (see Figure [2|). 
The consequence of a growing mode on super-horizon scales during inflation suggests that there is a deviation from 
scale invariance for large k. As a consequence, any observed deviation from scale invariance on these scales could be 
used to constrain these models. 

If we consider a Friedmann regime (—4.73 < r < —3.69, —1.27 < r < —0.81), the perturbations on large scale grow 
as if the effective fluid was a form of standard matter. For example, supposing that the expansion of the background 
follows the standard behaviour a cx i^/s^ which corresponds to r « —4.17, the parameter p(r) of the scalar field is 
« 3.6, and in this case w w 0. If we choose r = —4, we would have a radiation-like effective fluid etc. This is means 
that a non-minimally coupled scalar field in a Friedmann regime can mimic baryonic fluids, but do not necessarily 
interact with photons via Thompson scattering. 

Another curious behaviour in the Friedmann regime is that, for (—4.73 < r < -4.17,-1.07 < r < —1.27), the 
perturbations are effectively frozen. This feature is very important when one makes the transition from perturbations 
seeded in the early universe to the matter power spectrum (if we take (p to correspond to the inflaton) , and in general 
will affect the spectrum of scalar perturbations in the matter dominated era. 
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2 4 6 8 10 

FIG. 1: Plot of the solutions (|85|l in the inflationary regime (case r — 0.01). The solid curve represents the 'I'/^'o the dashed 
one represents A/^Pq. 




2 4 6 8 10 

FIG. 2; Plot of the solutions lf85|l for the case r = —5. The solid curve represents the 'i/^o the dashed one represents A/5'o. 



Finally, from the above equation we can also obtain the behaviour of the Newtonian potential ^n- 

{3 + rf (6 + 5r) 



N 



18 (12 + r (6 + r)) (6 + r (9 + 2 r)) 



C 



2 (3 + r) (6 + 5r) 

3 6 + r (9 + 2 r) 



(87) 



As in the minimally coupled case, <i> has a constant mode, but unlike the minimally coupled case, the other mode is 

12 + 6r + 

not necessarily decreasing. In fact, the exponent r — is negative only if —3 < r < (see Figure This 

r(3 + r) 

means that, when r is outside the above interval, the gravitational potential increases in time so that the rate at 
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FIG. 3: Plot of the exponent of the non-constant mode in (|87[l 



which the perturbations grow is faster than in the case for a minimally coupled scalar field. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have developed the theory of cosmological perturbations of a FLRW background for a generic 
scalar-tensor theory of gravity in vacuum. Our analysis is based on the fact that the field equations of a scalar-tensor 
theory in vacuum can be recast in a form that resembles the standard Einstein theory, in which an effective fluid is 
present. This effective fluid is perfect in the background (like in the case of a minimally coupled scalar field), but it 
is imperfect in the perturbed universe. This is due to the peculiar form of the heat flux and the anisotropic pressure 
that depend on the acceleration and the shear as a result of the non-minimal coupling F . 

As in pT| . given our choice of u"" , the perturbations can be characterized using the variables Da = ji^^aV afJ^, 
Ca ~ a?V aR and — ip~^aWa'ip- The first one represents the density perturbation of the effective fluid; the second 
one, the projected gradient of the 3-curvature scalar of the (j) — constant surfaces respectively; and the third one 
represents the inhomogeneities of the scalar field. Note that, due to the fact that in our frame Vq0 — 0, this variable 
depends on the momentum ■(/;. 

The three variables above are connected to each other via the constraint (|42p . In the minimally coupled case such 
constraint reduce to a proportionality relation between Da and \E'a and is used to eliminate the variable '^a- However, 
in our case ^'^ and Da have a rather different meaning and different relevance depending on the problem one wishes 
to solve. There can be cases in which the the perturbations of the effective fluid have a different impact on the matter 
perturbation and it is more useful to deal with Va than ^P^. 

Once we have chosen the perturbation variables it is relatively easy (although rather long!) to derive their evolution 
equations. These equation can be simplified if we consider spacetimes that do not differ too much from the background. 
In fact, in this case one can use the a linearization procedure based on the structure of the background to completely 
drop the interaction terms. 

In this paper, we limited ourselves to analyze spherically symmetric clumping. This has been done by considering 
only the projected divergence of the variables above. The equations for these variables describe the evolution of scalar 
perturbations in a general scalar-tensor theory of gravity in vacuum. They are valid also for any spatial geometry and 
can be used to analyse the perturbations in any scalar - field - dominated universe models (not necessarily inflationary). 

As an example we used a coupling, a potential and a background solution obtained in the framework of the 
Noether Symmetry approach (3^ . We were able to solve the perturbation equations exactly for this background 
in the long wavelength limit. Our solution reveals that there exist power law inflation regimes in which the scalar 
field perturbations grow instead of being dissipated. Such a behaviour will affect the matter power spectrum, and 
potentially break scale invariance. However, this conclusion depends strictly on the behaviour of the perturbations 
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on sub-horizon scales, which requires a complete analysis of the perturbation equations for all scales. 

Also, we found that in "Friedmann regimes" the effective fluid associated with a non-minimally coupled scalar field 
is able to mimic a standard perfect fluid that does not necessarily interact with photons. The perturbations of this 
fluid grow for most of the values of r associated to a Friedmann evolution, but there are also cases in which the 
super-horizon perturbations are effectively frozen, i.e., their growth rate is slower than the expansion rate, thereby 
differing significantly from the standard gravitational instability picture. 

Finally, we derived the evolution of the gauge invariant version of the Newtonian potential (/>jv. As in the case of 
a minimally couple scalar field, (pN contains a constant mode. However, unlike the minimally coupled case, it also 
admits a second mode that is not necessarily decaying and this again modifies the standard picture. 

In conclusion, our preliminary results reveal that the presence of a non-minimal coupling has unexpected and 
interesting effects on the evolution of the long wavelength scalar perturbations. An analysis of the full perturbation 
spectrum will clarify further on these properties and will be presented in a forthcoming paper. 
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For sake of simplicity in the main body of the paper we give here the expressions of the coefficient of the equations 
(j6ip . All the quantities appearing in these coefhcients are meant to be the background ones. 
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3F(l + w)M'/'^ + 0^(2V'^-3F(l + w)/i-2e^^ + 3F 



BA[t) 



KF^ 3 F/i (2^92-3 (1 + w) ^) - 6 6^ _^ 9 F (2 6^ -9^i)-QQ^F 



+ 




(Al) 



Ca(0 




2KeF^ 49F-3V'M2 + F") 




(A2) 



-DAit) 
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3F {1 + w) (2Fi;'^+3F'^'^ 

6F2 + iiFtl)'^ + QF@F'^ {2ip'^-F {1 + w) ji) 
3F-^ {1 + w) jj? 



F^ (qF {1 + w) iJ. + 12eF -^FqA +6FeF'^tp'^ F" 

+ ^ 7 ^—^ (A3) 

3F'^ {1 + w) f2FV2 + 3F2j 

_ 3FM(l + ^;) + 3e^ 

-^^^'^ - -l^F^^eliTwr ^ ^ 

6F II {2we'^ -3 {1 + w) ij) -ue^i^^ + 2eF {20^ -ue^ F 

^^^'^ = ZF{l+w)Qi, ^^^^ 



49 (SV'^ (2 + i?^") -49^) 



2a2F (292 + 3/x) +120^9 (^F w n - i/j'^ - P'^ 



2 a 



J^F (9 {1 + w) /i-892) +129F2 + 6F9V^ (2 + F") 



= 9,F2 (1 + .) ' 

gJsTF (f (1+w;) ii + &f) 
A^{t) = — ^ ' 



+ 



2a^ {1 + w) (2FV2 + 3F2^ 
6/x9F2 {l + w)+F [f(492 + 9 {1 + w) n) +69F 



- 6F9F 



B^{t) 



+ 



4 0^92 (2FV'2 + 3F2) 

4 // F- B (2 (cf- + () // + 9 /(■ /() - 18 iiFQF 

6 92 (^2FV'2 + 3F2j 

3/iF (2F92 + 9F (l + w) /i + 69VF') 

692 (^2FV'2 + 3F2^ 

F' (3Fy^ (2^92 -3 (1 + w) At) -692^2) 
a2 {1 + w) 92 /i^ (2F + 3F'2) 

2 V - 9 V) - 6 F' 92 F 
a2 {1 + w) e^n (2FV2+3F2) 



(AlO) 



^ 9F^ (292^-9mV')-6F^9^F ^^^^^ 
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2KF (^QF-Sij^ {2 + F") 



a? (l + w) /X (^2FV'2 + 3F2 

F^j^ [2 62-3 (/i + Sw^i)] +39 [aF ^j^' + ?, F^ Q 
36 (^2FV'2 ^_3/^2 

36i^ (2i^ (62+/i + 3?i;Ai) +3F 



36 



(2Fij^ 



3F2 



(A12) 



2^a26 (l + w) [2Fi:^ + 3F^ 



(A13) 



6 6 F (t/-^ _ ^ ^ ^ j^'^ _ ^2 (2 6^ + 3 /i) 



36 



(l + w) (2FV'2^3/^2j 



(A14) 



2i^ 



3F (F (l + w) ^-2V'2) + 4F6i^-6i^2_3^y/^// 



3/iF (l + w) (2FV'2 + 3F2 



(A15) 
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is not a conserved quantity for scales larger than Hubble radius [ll[. For this reason we will not include it in our discussion. 
[38] Note that the underlying assumption in this decomposition is that the perturbation variables can be factorised into purely 

temporal and purely spatial components. 
[39] This solution can represent, of course, also a contraction, but in the following we will consider only the expanding cases. 
[40] This is a peculiar feature of this background solution that is not true in in general. In fact it can be easily verified that it is 

originated by the peculiar relation between the time behaviour of the scalar field and the power appearing in the potential 

of this example. 



